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Abstract 

In this work, we study the 'scalar channel' of the emission of Hawking radiation from a 
(4+n)-dimensional, rotating black hole on the brane. We numerically solve both the radial 
and angular part of the equation of motion for the scalar field, and determine the exact 
values of the absorption probability and of the spheroidal harmonics, respectively. With 
these, we calculate the particle, energy and angular momentum emission rates, as well 
as the angular variation in the flux and power spectra - a distinctive feature of emission 
during the spin-down phase of the life of the produced black hole. Our analysis is free 
from any approximations, with our results being valid for arbitrarily large values of the 
energy of the emitted particle, angular momentum of the black hole and dimensionality of 
spacetime. We finally compute the total emissivities for the number of particles, energy 
and angular momentum and compare their relative behaviour for different values of the 
parameters of the theory. 



1 Introduction 



The formulation of the theory with Large Extra Dimensions P was driven by the motiva- 
tion to address the hierarchy problem (for some early works, see j2]). While all Standard 
Model fields are restricted to live on a (3+l)-dimensional brane, gravity is allowed to 
propagate in the (4 + n) -dimensional bulk; above the scale of decompactification of extra 
dimensions M*, that can be as low as 1 TeV, gravitational interactions become strong 
and thus comparable to the Electroweak interactions. Moreover, above M,,, all effective 
theories cease to be valid, and a more fundamental theory, that describes all forces in- 
cluding gravity, must take over. In the context of such a theory, transplanckian collisions 
of particles become particularly important as their product cannot be ordinary particles 
any more but rather heavy objects. This opens up the exciting possibility of observing 
strong gravitational phenomena, linked to a Quantum Theory of Gravity, at a low energy 
scale, possibly at the TeV scale. 

The creation of mini black holes as the result of such transplanckian particle col- 
lisions [S] has drawn considerable interest during the last few years. Such black holes 
may be created either at colliders jl] or in high energy cosmic-ray interactions [HI (for 
an extensive discussion of the phenomenological implications and a more complete list of 
references, see the reviews [HI El IE])- If the produced black holes have a mass consider- 
ably larger than the fundamental Planck mass M^,, quantum gravity effects can be safely 
ignored, and the black holes can be treated as classical objects. Their most important 
characteristic, as well as the most prominent signature of their creation, will be the emis- 
sion of Hawking radiation in the form of elementary particles. This will take place both 
in the bulk and on the brane, with the latter 'channel' being the most important from the 
phenomenological point of view. After shedding all additional quantum numbers during 
a short balding phase, as dictated by the no-hair theorem of General Relativity, the pro- 
duced black hole will settle down to a Kerr-like phase, the spin- down phase, during which 
the black hole will mainly lose its angular momentum through the emission of Hawking 
radiation and superradiance. After that, the Schwarzschild phase will commence with the 
black hole (now spherically symmetric) emitting Hawking radiation and gradually losing 
its actual mass. 

The Schwarzschild phase of the life of a small higher-dimensional black hole formed 
in a flat background has been admittedly the most well studied in the literature. The 
task of determining the spectrum of Hawking radiation during this phase has been at- 
tacked both analytically jUJ El and numerically ^2]- While the first set of works 
led to the derivation of useful, analytical formulae for the emission rate, the latter work 
provided exact, numerical results valid at all energy regimes. Both approaches reached 
the conclusion that the Hawking radiation spectrum strongly depends on the existence of 
additional spacelike dimensions in nature, with this dependence being reflected both in 
the number of particles and energy emitted by the black hole per unit time as well as in the 
actual type of particles emitted (scalars vs. fermions vs. gauge bosons). Further studies 
have also shed light on the dependence of Hawking radiation from a higher- dimensional 
black hole on higher-derivative gravitational (Gauss-Bonnet) terms 13^, the mass of the 
emitted particles OUni, the charge of the black hole [16], and last but not least, the 
cosmological constant [T7], with the latter leading to a distinct signature of its existence 
at the low-energy part of the radiation spectrum. 
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The progress in studying the emission of Hawking radiation during the spin-down 
phase of a small, higher- dimensional black hole has been much slower. This was due to the 
more complicated gravitational background but also to additional technical difficulties in 
solving the equation of motion for a particle propagating in such a background. Analyti- 
cal formulae for the emission rate were derived flEl IT^ but the results were partial, being 
valid only at the low-energy regime, for low angular momentum of the black hole, and for 
a specific dimensionality of spacetime. The first exact numerical results for the Hawking 
radiation emission rate in the form of scalar fields, during the spin-down phase, as well 
as for the amplification due to superradiance, were presented in [201 • These were valid 
for arbitrary values of the energy of the emitted particle and for arbitrary values of the 
angular momentum parameter of the black hole. That was the result of the determination 
of the exact value of the angular eigenvalue - being also involved in the radial part of 
the equation of motion of the particle - via numerical means instead of the use of an 
approximate analytical formula. The analysis in demonstrated that the suppression 
of the energy emission rate at the low-energy regime as the angular momentum of the 
black hole increases, as was found in is in fact overturned in the intermediate- and 
higher-energy regime where a strong enhancement appears instead. In the case of superra- 
diance, these results confirmed the existence of the effect in the higher-dimensional case, 
and demonstrated the exact energy amplification: although strongly mode-dependent, 
this amplification can reach approximately a 10% magnitude^. Shortly afterwards, an- 
other study appeared in the literature focused on the study of the power (energy) 
spectrum, as well as its angular distribution in the 5-dimensional case, but once more the 
analysis relied on the assumption of low energy and low angular momentum. Two more 
studies focused on 5-dimensional rotating black holes appeared recently [211 ESj, where 
the question of the emission in the bulk was also addressed. 

At present, a comprehensive study of the emission of Hawking radiation from a 
higher-dimensional black hole in its spin-down phase, that could provide exact numerical 
results for arbitrary values of the angular momentum of the black hole, the energy of the 
emitted particle and the dimensionality of spacetime, is still missing from the literature. 
Apart from the power (energy) spectrum, whose study needs to be completed for arbitrary 
values of the aforementioned parameters, the question of the emission rate of particles 
(flux spectrum) and of the rate of loss of angular momentum of the rotating black hole 
also needs to be addressed. Moreover, the angular distribution of the emitted radiation - 
a distinctive signature of emission from a rotating higher- dimensional black hole - needs 
to be investigated, with the results again being free from any restrictive assumption about 
the values of the fundamental parameters of the theory. It is these tasks that we have 
undertaken to fulfil in this work. The strongly spin-dependent techniques needed to be 
apphed for the calculation of the spectra for different types of particles, and the plethora of 
results that need to be derived in each case, have forced us to restrict our attention in this 
work to the case of scalar fields emitted by the black hole, and present the corresponding 
analysis and results for higher-spin fields in a subsequent work [22] • 

Before presenting the outline of our paper, a brief discussion of the assumptions 
made during our analysis should be added here. As mentioned earlier, it will be assumed 
that the black hole mass Mbh is considerably larger than the fundamental scale of gravity 
in order to ensure that quantum effects are small and that the black hole geometry 

^Results on the superradiance effect in the case of a higher-dimensional black hole were also presented 
in EI] and [221- 
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can be safely considered as a classical object. The same constraint guarantees the absence 
of any significant back reaction to the gravitational background due to the change in the 
black hole mass after the emission of a particle: the assumption that the temperature of 
the black hole, around which the emission spectrum is centered, is much smaller than the 
black hole mass translates again to Mbh ^ We are also going to assume that the 
horizon of the black hole r/j is much smaller than the (common) size L of the additional 
compact spacelike dimensions so that the black hole can be considered as embedded in 
a (4 + n) -dimensional, non compact, empty spacetime. Finally, since the value of the 
brane self-energy can be naturally assumed to be of the order of the fundamental Planck 
scale M*, and thus much smaller than the black hole mass, its effect on the gravitational 
background can be ignored. 

We will start, in Section 2, by presenting the theoretical framework for our analysis 
including basic formulae for the gravitational background, the equation of motion for the 
propagation of a scalar field in it, and the Hawking radiation emission spectra. In section 
3, a description of the numerical methods and techniques used for the determination 
of the exact value of the angular eigenvalue and of the angular and radial part of the 
wavefunction of the field will be given. After the angular eigenvalue is determined, the 
radial part of the equation of motion is numerically solved to determine the absorption 
probability (greybody factor) and thus the various spectra. In section 4, we present our 
results for the flux, power and angular momentum spectra (integrated over the azimuthal 
angle 6) for arbitrary values of the energy of the particle, the angular momentum of the 
black hole, and the dimensionality of spacetime. All the spectra have a thermal profile in 
terms of the energy of the particle, as expected, and bear a strong dependence on the value 
of the angular momentum of the black hole and the number of extra dimensions. Next, 
we numerically solve the angular part of the equation of motion to find the exact values 
of the spheroidal harmonics, and we calculate the angular distribution of the flux and 
power spectrum; these results are also given in section 4, and, contrary to the emission of 
radiation during the Schwarzschild phase, are characterized by a strong angular variation. 
Finally, in the last part of section 4 we focus on the study of the various total emissivities 
(the total number of particles, energy and angular momentum emitted per unit time by 
the black hole in all frequencies) and compare their behaviour for different values of the 
dimensionality of spacetime and of the angular momentum of the black hole. We finish 
with a summary of our results and conclusions, in Section 5. 



2 Theoretical Framework 

The gravitational background around a (4 + r;,)-dimensional, rotating, uncharged black 
hole was found by Myers & Perry [27j, and in what follows we will use this line-element 
to describe the spin-down phase of a small, higher- dimensional black hole created by 
the collision of highly energetic particles. If we assume that the colliding particles are 
restricted to propagate on an infinitely-thin 3-brane, they will have a non-zero impact 
parameter only along our brane, and thus acquire only one non-zero angular momentum 
parameter about an axis in the brane. Such a higher-dimensional black hole will then be 



where 



described by the following line-element [2Z] 

+ + sin2 Q dip^ - cos^ 6 dQn, (1) 

A = r2 + a2--^, J: = r^ + a^ cos^e, (2) 

and is the line-element on a unit n-sphere. The mass and angular momentum (trans- 
verse to the ry9-plane) of the black hole are then given by 

{n + 2)A„+2 , 2 

Ms/^ = — — /X , J = — — - Mbh a , (3) 

IottG n + 2 

with G being the (4 + n)-dimensional Newton's constant, and A„_|_2 the area of a (n + 2)- 
dimensional unit sphere given by 

27r("+3)/2 

= r[(n + 3)/2] • 

In this work, we will focus on the emission of Hawking radiation [2H], in the form of 
scalar fields, from this (4 -|- n)-dimensional, rotating black hole ^. For phenomeno logical 
reasons, we will concentrate on the emission of these fields directly on our brane since any 
particle modes emitted in the bulk cannot be detected by an observer restricted to live 
on the brane. We will thus need to determine first the line-element on which the brane- 
localized modes propagate, and then to solve their equation of motion on the resulting 
background. The induced-on-the-brane line-element can be found by fixing the values of 
the additional angular coordinates that were introduced to describe the compact extra 
n dimensions: by setting 6i = 7i/2, for i = 2, ...,n + 1, we are led to the 4-dimensional 
background 

ds^ =(l- -J^) de + dt d^ - ^dr^ 

-Ede^-(r^ + a^ + ^^^^) sin^^V. 

The black hole horizon is given by solving the equation A(r) = 0, which, for n >1, 
leads to a unique solution given by 

where we have defined a* = a/vh- In addition, while for n = and n = 1, there is 
a maximum possible value of a, that guarantees the existence of a real solution to the 
equation A = and thus of a horizon, for n > 1 there is no fundamental upper bound 
on a and a horizon always exists. An upper bound can nevertheless be imposed on 



^For some classic works on Hawking radiation in the 4-dimcnsional spacetime, see |29l l3(Jll^ . 
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the angular momentum parameter of the black hole by demanding the creation of the 
black hole itself from the collision of the two particles. The maximum value of the impact 
parameter between the two particles that can lead to the creation of a black hole was 
found to be El 



1 + 



. 2- 1 



, (7) 



an analytic expression that is in very good agreement with the numerical results produced 
in jn]c. Then, by writing J = hMBn/'^i for the angular momentum of the black hole, and 
using Eq. ^ and the second of Eq. Q, we obtain 



^max _ n + 2 



The equation of motion for a particle with spin 0, 1/2 or 1, propagating in the 
induced-on-the-brane gravitational background (0) of a higher-dimensional, rotating black 
hole was derived in jHlEj. For scalar fields, the field factorization 

r, e, if) = e-^"* e*'"^ R{r) Tp{e, au) , (9) 

where Tp{9,auj) are the so-called spheroidal harmonics jHSj, leads to the following set of 
decoupled radial and angular equations 



1 d ^ drne^ ^ + + cos' 9 + E^^ Trie, aco) = 0, (11) 



sin 9 d9 

respectively. In the above, we have defined 

K = {r' + a') uo-am, = + a'uj'^ - 2amuo . (12) 

The angular eigenvalue E'^ provides a link between the angular and radial equation, and 
various methods for finding its exact value will be discussed in the next section. Both 
equations (fTm) -(fTT j) will be solved numerically to determine the radial and angular parts 
of the scalar field. The radial equation (|l(jp can be solved analytically in the asymptotic 
regimes of near horizon and infinity. By using the new radial, 'tortoise', coordinate 

^ = (13) 

dr A(r) 

we find the following asymptotic solution near the horizon of the black hole 

i?^(r*) = e'^'' + A2 e-'^'* , (14) 



where 2 are integration constants, and k is defined as 



A boundary condition must be applied in the near-horizon regime to ensure that the 
solution contains only incoming modes; this is satisfied if we set Ai = 0. On the other 
hand, for fixed a* and large r, the asymptotic solution at infinity takes the form 

R^{r) = B, — + B2^, (16) 

where -81,2 are again integration constants. The asymptotic solutions (fT^ and (fTHjl will 
serve as boundary conditions for our numerical analysis. 

The Hawking temperature of the (4 + n)-dimensional, rotating black hole is found 

to be 

_ (n+l) + (n-l)a^ 
^» - Ml + al)r, ' 
and leads to the emission of thermal Hawking radiation both in the bulk and on the brane. 
As mentioned earlier, here we will study the scalar field 'channel' emitted on the brane. 
Since we are interested in studying an evaporating black hole, the relevant quantum state 
is the "past" Unruh vacuum \U~) |3Sl OH tHHI- This state is defined in terms of the 
standard "in" and "up" modes, corresponding to unit initial particle fiux at past null 
infinity and the past horizon, respectively [3S1IM1ISSI- At infinity, far from the black hole, 
the outward fiuxes of energy E and angular momentum J can be expressed in terms of 
the components T*"* and of the renormalized stress-energy tensor for a quantum scalar 
field in the state |f/-) jSll: 

dE f 

{U-\T'^'\U-)rcnr\med9dip, 



dt 
dJ 
lit 



{U-\T;^\U-),er^r^ sine d9d^, (18) 
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where the integral is taken over the sphere at infinity. These two components do not 
need renormalization for a quantum scalar field on the metric (0) [34 , thus simplifying 
the computation. Our boundary conditions, (fT^ and (fT^ . correspond to the "in" modes 
[33. , and in terms of these the relevant quantities can be written as^ 

iu-\Tn. - i^E£„557|^;|pTM-nTO-)r. (19) 

<P-\T^\U-U - E £ S^^i^^lAJ^ |Tr(9.o.)P . (20) 

In the above expressions, |^^,mP is the absorption (or transmission) probability for a scalar 
particle propagating in the brane background (0), and k is given by ()15|) . Multiplying 
these expressions by sin 6 and integrating over the sphere at infinity gives the total 
outgoing fiuxes of energy and angular momentum: 

~Tr = 7^52 [ r, /rj. 1 7|^<?,mP, (21) 

dt 27T y^^o exp [k/Tu] - 1 



dJ 1 ^ r°° mduj 



dt 27r f - £=0 exp [/c/Th 



^l,m\ 1 (22) 



■^Wronskian relations between the "in" and "up" modes allow us to easily change our basis; in terms 
of the "up" modes, the quantity \At,m\'^ = l-A™^^ in Eqs. I^Hll-lEni is replaced by {u)/k) \B~^^f in the 
notation of \6'6\. 



where we have used the fact that the spheroidal harmonics T^{6, auj) are normahzed 
according to 

/•TT 

de sme\T^{e,auj)\^ = 1. (23) 
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The absorption probabihty |.4.^,m| has an exphcit dependence not only on the angular 
momentum numbers {i,m), as denoted, but also on the energy u of the emitted particle 
and the number of extra dimensions n. Therefore, its presence in the expressions of 
the emission rates modifies significantly the blackbody profile of the spectrum, especially 
in the low- and intermediate-energy regime. Its exact form will be found by solving 
numerically the radial equation (llUj) and determining the amplitudes of the incoming 
and outgoing modes at infinity, according to the asymptotic solution |T6|l . Then, the 
absorption probability may be written as 



5l 



(24) 



The first part of this work will focus on the computation of the differential emis- 
sion rates of particles (flux spectrum), energy (power spectrum) and angular momentum, 
integrated over all angles Q. These are given by 

d^N 1 1 



dtduj 2tx ^ exp [/c/Th 
d^E 1 V- uj 



dtduj 27r ^ exp [A;/Th 
£■.] 1 \ - m 



dtdu) 2tt ^ exp \klT-a\ — 1 



By using the exact solutions of the angular equation (fTTj). we will then study the angular 
distribution of the differential fluxes of particles and energy due to the axial symmetry 
of the gravitational background - a feature that will be a distinct signature of emission 
during the spin-down phase: 



d^N 



d{cos9)dtdLU 27r exp [/c/Th] — 1 



d^E 



d{cos6)dtduj 27r exp [/c/Th 

In both cases, the various spectra will be given as a function of the energy of the emitted 
mode u, but their dependence on the angular momentum parameter of the black hole and 
the dimensionality of spacetime will be also studied in detail. 



3 Numerical Analysis 

For the purpose of the analysis presented in this work, we will need to solve both the 
angular (lllj) and the radial ()10j) part of the scalar equation. The angular equation will be 
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studied first as the eigenvalue E'^ is required before we can integrate tlie radial equation. 
The spheroidal harmonics, which are the solutions of Eq. have been extensively 

studied in the literature [32 , therefore, here we will only briefly outline our method and 
some of their key properties. 

Before describing the numerical analysis performed in this work though, we would 
like to note that if one is interested in solving the radial equation only and studying the 
spectrum integrated over the angle 6, then it is sufficient to find only the eigenvalues E^, 
while the eigenfunctions T^{6, auj) are not required. The eigenvalues E^ are functions of 
auj, and an analytic form exists only in the limit of small auj. To derive results for the 
emission of Hawking radiation valid for arbitrarily large energy of the emitted particle 
and angular momentum of the black hole, one needs to find numerically the exact value 
of EJ^. This can be done, for instance, by using the so-called continuation method |36j . 
The continuation method is a generalization of perturbation theory which is applicable 
when the changes in the initial Hamiltonian (for which the eigenvalues are known) are not 
necessarily small. It involves writing the T^{6, au) functions in the basis of the 6'-parts 
of the spherical harmonics, SY^{6), which in conjunction with the angular equation and 
techniques from the perturbation theory leads to a differential equation for the eigenvalue 
E^. This is, then, solved numerically, by using appropriate initial conditions, to derive 
the eigenvalue for any i and m and for any value of au. This method was outlined in 
jSj and used in pUj for the derivation of the exact power spectrum for the emission of 
Hawking radiation, in the form of scalar fields, from a 5-dimensional black hole on the 
brane. 

As mentioned above, in this work we also compute the exact angular distributions 
of the fiuxes of particles and energy, therefore we require the values of the spheroidal 
harmonics T[^{6,auj) themselves. For this, we follow a modified version of the method 
presented in section 17.4 of [37J. Given the symmetry in the angular equation (lllj) under 
the replacement of m by — m, it is sufficient to consider only those solutions for which 
m > 0. It is also convenient to change the independent variable to rj = cos 6 and define a 
new dependent variable yimuj by 



(1 - V') -f^ - 2 (m + 1) + [ET - m (m + 1) + a^uj^^) y,m. = . (31) 



The above equation must be solved on the interval [—1, 1], with the boundary conditions 
that yemuj is regular and non-zero at both end-points. We used a shooting method jHZ| to 
solve this boundary value problem. For spin particles, the number of zeros of yimuj in 
the interval [—1, 1] is given hj i — m, and the spheroidal harmonics are either odd or even 
functions of rj. This symmetry can be used to integrate Eq. 1)311) over the half-interval 
[— 1, 0], using appropriate boundary conditions at ?7 = 0. However, we have chosen not to 
take this path, but instead to integrate over the full interval [—1,1] as this method will 
more readily generalize to higher spins [211 , when the spin- weighted spheroidal harmonics 
do not have this symmetry property. 

The end-points, i] = ±1, are regular singular points of Eq. (jHH). By using the 
Frobenius method, we find two linearly independent solutions near these end-points, one 



ytrnu) 



(30) 



in terms of which the angular equation (|TT|) becomes 
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Figure 1: Spheroidal harmonics for a scalar field with i = m = 0, and for auj = (0, 1, 2, 3). 



of which is regular and the other irregular. We begin the numerical integration at ?7 = 
— 1 + 61], with 6t] = 10~^, using the Frobenius series expansion to impose our initial 
condition for a regular function y. We then integrate out to 1] = 1, and demand that the 
solution be also regular at 17 = 1. For a general initial value of E^, the numerical solution, 
close to ?7 = 1, will be a linear combination of the regular solution we are seeking and the 
irregular solution. We define a function F, which depends on E^, as follows: 



where y is the numerical solution we have found by integration, and j/reg is the solution 
regular at 77 = 1. The idea is therefore to find the value of E^ such that F{E^) = 0. Due 
to the irregular solution, the integration cannot be extended all the way to r/ = 1, but 
instead we impose the boundary condition at = 1 — 61], again using a Frobenius series 
expansion for the regular solution y^^g- 

The above method can give us the values of both the angular eigenvalues E^^ and 
the spheroidal harmonics Tf^{6,auj). In the case of the eigenvalues, we have used both 
the continuation and the shooting method, and we have found an excellent agreement 
between the values obtained. For the spheroidal eigenf unctions, the results found by 
using the shooting method^ were finally normalized according to Eq. (j23|) . In Fig. Q as 
an illustrative example, we depict the derived spheroidal harmonics for the scalar mode 
with £ = m = 0, and for the values aoo = 0, 1,2 and 3 - we remind the reader that, 
for au = 0, the spheroidal harmonics reduce to the spherical ones. In previous works 
[T^ I2S] , the spherical harmonics were used as an approximation to the exact spheroidal 
ones. From Fig. Q we can clearly see that, as long as atu <^ 1, this is indeed a valid 
approximation. However, as au increases, this approximation becomes increasingly poor. 
In addition, the difference between the spheroidal harmonics and the spherical ones is 
strongly mode- dependent, being more significant for the low {i,m) modes and less so for 
the high {i, m) modes. For an accurate analysis therefore, leading to the exact angular 

^Further details of this numerical method, generalized for all spins, can be found in 23 ~ an alternative 
method of finding the spheroidal harmonics, once the eigenvalues are known, is presented in |39| . 
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distribution of particle and energy fluxes, valid at all energy and angular momentum 
regimes, the use of the exact spheroidal harmonics is imperative. 

Having derived the values of the angular eigenvalues E^, the integration of the radial 
differential equation (fTIH) can now proceed. We start the integration at the horizon of the 
black hole and integrate outwards towards infinity. Comparing our numerical results 
for the radial solution with the asymptotic solution at infinity (fTB|) . we determine the 
integration constants Bi and B2. The absorption probability for scalar fields can then be 
derived from the relation (j23). The use of the exact value of the angular eigenvalue E^, 
instead of an approximate analytical formula, guarantees the validity of the solution of the 
radial equation, too, for arbitrary energy of the emitted particle and angular momentum 
of the black hole. 

4 Numerical Results 

In this section, we present exact numerical results for the flux, power and angular momen- 
tum spectra for the emission of scalar particles on the brane from a higher-dimensional, 
rotating black hole. The dependence of the various spectra on the two fundamental pa- 
rameters - the dimensionality of spacetime n and the angular momentum parameter a^, - 
will be studied in detail, and in the case of the first two types of spectra (flux and power), 
the angular distribution of the emitted particles and energy will also be examined. We will 
finally study the corresponding total emissivities of the black hole on the brane. These 
follow by integrating the various emission spectra over the whole frequency regime, and 
can provide information on the dependence of the total number of particles, energy and 
angular momentum, emitted by the black hole per unit time, on n and a^, as well as on 
their relative behaviour as the values of n and a^, vary. 

4.1 Flux emission spectra 

We start by presenting our numerical results for the number of scalar particles emitted 
by the black hole on the brane per unit time and unit frequency. Figures I2ta,b) depict 
the flux emission rate, Eq. ()25|) . as a function of the frequency of the emitted particle 
in dimensionless units^, oovh, for various values of the angular momentum parameter 
of the black hole. The behaviour of the flux depends strongly also on the dimensionality 
of spacetime, and, for that reason. Figs. |2fa) and (b) present results for two indicative 
cases, n = 1 and n = 4, respectively. 

In the case of a 5-dimensional black hole, the particle flux, shown in Fig. |2ta), is char- 
acterized, for small values of a,,, by a strong peak at low values of uvh, this clearly favours 
the emission of low-energy quanta. As the angular momentum parameter increases, how- 
ever, this peaked Gaussian curve gives its place to a broad one, where high-energy particles 
become almost as equally likely to be emitted as the low- and intermediate-energy ones. 
This feature will be of prime importance in the determination of the power spectrum, to 
be studied in the next subsection. For higher values of a^,, oscillations become a com- 
mon feature of the spectrum; these are due to the higher partial waves gradually coming 

^Throughout this analysis, we will be assuming that the horizon radius rji of the black hole remains 
fixed as n and a, varies. 
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Figure 2: Flux emission spectra for scalar particles on the brane from a rotating black 
hole, for (a) n = 1, and (b) n = 4, and various values of a*. 

into dominance as the energy increases. In Fig. Efb), we depict the flux spectrum of 
an 8-dimensional rotating black hole. The spectrum shares common characteristics with 
the one obtained in the 5-dimensional case, with the curve becoming broader and the tail 
dying away much more slowly, as a^, increases. However, there is now a clear enhancement 
in the number of particles emitted by the black hole per unit time over the whole energy 
regime, as a^, increases, and the oscillations, although still present, are less apparent here. 

Comparing the vertical axes of Figs. Efa) and (b), one easily observes the almost 
one order-of-magnitude enhancement in the flux of particles, as we go from the n = 1 to 
the n = 4 case. This enhancement was a characteristic feature of the flux spectrum in the 
case of a non-rotating black hole ^2]; and - as we see in Fig. it characterizes also the 
spectrum of a rotating black hole. Figure El depicts the particle emission rate for a rotating 
black hole with a fixed angular momentum parameter, = 1, living in spacetimes with 
different dimensionalities. The particle flux broadly increases, and peaks for higher values 
of ujTh, as n increases, clearly leading to a higher total number of particles emitted per 
unit time, as we will also see in subsection 3.4. 




Figure 3: Flux emission spectra for scalar particles on the brane from a rotating black 
hole, with = 1, for various values of n. 



11 



Flux - Flux Flux 




Figure 4: Angular distribution of the flux spectra for scalar emission on the brane from a 
rotating black hole, for n = 1 and = (0, 0.6, 1). 




Figure 5: Angular distribution of the flux spectra for scalar emission on the brane from a 
rotating black hole, for n = 4 and a^, = (0, 0.6, 1). 



Also of great importance is the angular distribution of the emitted particles. Where- 
as, in the case of a non-rotating black hole, we expect a uniform distribution of particles, 
in the case of a rotating black hole, we anticipate a distinctly different angular distribution 
pattern. An angular variation in the flux and power spectra would be a clear, characteristic 
signature for emission during the spin-down phase in the life of a black hole. Figures 0] 
and El depict the flux emission spectra on the brane from a rotating black hole with n = 1 
and n = 4, respectively, and for a* = (0, 0.6, 1), as a function of the energy of the emitted 
particle and the value of cos^ of the azimuthal angle. As expected, for vanishing a^,, the 
spectrum shows no angular variation, independently of the value of n. As the angular 
momentum parameter increases, the emitted modes are starting to concentrate on a region 
around the equator {6 = vr/2), however, their angular distribution depends, at the same 
time, also on their energy: even for non- vanishing a*, at low energies the spectrum is 
dominated by a spherically-symmetric distribution, and shows no dependence on the 
angle 6; only when the energy exceeds a certain value does the non-spherically-symmetric 
nature of the emission become significant. In general, for fixed a*, as n increases, the 
non-spherically-symmetric behaviour becomes dominant sooner, i.e. at lower energies, 
and there is an increasing contribution from angles away from the equator ^. 

^We should note here that after the end of our analysis, the angular distribution of both the particle 
and energy flux was integrated over cos 9 to check that this produces the original particle flux and power 
spectrum. The two results agreed to an accuracy of 10~^. 
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4.2 Power emission spectra 

We now move to the power spectrum that describes the energy emitted by the black hole 
on the brane in the form of scalar particles per unit time and unit frequency. As in the 
previous subsection, in Figs. inta,b), we display the energy emission spectra, for n = 1 
and n = 4 respectively, as a function of the frequency of the emitted particle and for 
various values of the angular momentum parameter a^, of the black hole. 

Since the energy emission rate follows by multiplying the particle emission rate by 
the frequency of the emitted particle, the numerical results shown in Figs. inta,b) can be 
easily justified by using the ones for the flux spectrum derived in the previous subsection. 
In the case of a 5-dimensional black hole, we found that, for small values of a^,, the black 
hole prefers to emit low-energy quanta, while, for higher values of a^,, the emission of 
low-energy particles is suppressed while that of intermediate- and high-energy particles 
becomes equally likely. Therefore, as the angular momentum parameter increases, we 
expect the energy emitted by the black hole per unit time to be suppressed in the low- 
energy regime, and to significantly increase in the high-energy regime - this is exactly the 
behaviour depicted in Fig. E^a). 

In the case of an 8-dimensional black hole, however, the flux spectrum [see Fig. I2tb)] 
shows no suppression even at the low-energy regime. As a result, the energy emission rate 
is uniformly enhanced, as a^, increases, over the whole frequency regime. This behaviour, 
shown in Fig. IHl^b), seems to be in disagreement with the one derived in Ref. [22]; where 
the suppression at the low-energy regime in the power spectrum persists for all values - 
low or high - of the dimensionality of spacetime. 

Other features characterizing the flux spectra are also transferred almost intact to 
the power spectra. The oscillations make again their appearance, as the angular momen- 
tum increases, being quite prominent in the n = 1 case and smoothed out for n = 4. 
In both cases, the emission curves peak at higher values of for larger values of a^^, 
signifying a clear enhancement of the total emission rate, as we will see in detail in Sec- 
tion 3.4. By comparing again the vertical axes of Figs. Efa) and (b), an enhancement of 
more than one order-of-magnitude in the energy emission rate of the black hole is easily 
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Figure 7: Power emission spectra for scalar particles on the brane from a rotating black 
hole, with a* = 1, for various values of n. 



observed, as we go from n = 1 to n = 4. The enhancement in the power spectrum, as 
the dimensionality of spacetime increases, is clearly associated with the one found for the 
flux spectrum in the previous subsection, and is shown in detail in Fig. [3 for fixed a* and 
various values of n. A similar, strong enhancement in the power spectrum was also found 
in the case of a non- rotating black hole in ^21, and independently confirmed by other 
subsequent works [T51 ^21 ^] ■ In the case of higher- dimensional, rotating black holes, the 
only other available results in the literature are the ones presented in however, the 
figures in that paper show only an infinitesimal enhancement of the power spectrum as n 
increases, and curiously enough the enhancement is also missing from the results in 
for the non-rotating, limiting case with a^, = 0. 

In Figs. IHl and IHl we depict the angular distribution of the power emission spec- 
tra from a rotating black hole on the brane for n = 1 and n = 4, respectively, and 
for a,, = (0,0.6,1). Again, for vanishing a*, the spectrum shows no angular variation, 
but, as the angular momentum parameter increases, the emitted modes are starting to 
concentrate on the region around the equator. At low-energies, the spectrum is again 
dominated by spherically symmetric behaviour, and thus remains independent of the an- 
gle 6, but gradually the non-spherically-symmetric emission becomes important and an 
angular variation appears in the spectrum. As a^, increases, we can clearly observe the 
enhancement in the energy emission rate as well as the shift of the peak of the curve 
towards higher values of utrh- Similar behaviour, although on a much bigger scale, can be 
seen as we go from the n = 1 to the n = 4 case. 

The angular distribution of the power spectrum of a 5-dimensional black hole (n = 1) 
was also studied in Ref. [23] ■ As mentioned before, these results were only approximate, 
in the sense that the spherical harmonics were used instead of the exact spheroidal ones 
(which is valid only for acu -C 1). It would therefore be useful to compare those results 
with the ones derived here, and thus check the range of validity of the approximations 
made in Ref. [SB]- By comparing our Figs. IHta,b,c) (up to urh = 1) with Figs. 8(a,b,c) in 

we may see that an agreement arises in general terms: as a* increases, the emitted 
modes start to concentrate near the equator, oscillations appear, and the peak of the 
emission curve moves toward the right. Some differences though do appear, too: (i) in 
our case, the spherically symmetric emission becomes sub-dominant much faster, and 
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Figure 8: Angular distribution of the power spectra for scalar emission from rotating 
black holes, for n = 1 and a* = (0, 0.6, 1). 




Figure 9: Angular distribution of the power spectra for scalar emission from rotating 
black holes, for n = 4 and a* = (0, 0.6, 1). 



the low-energy peak is rapidly suppressed, in contradiction with the behaviour depicted 
in Fig. 8 of Ref. [22]; (h) the height of the peaks appearing at higher energies are also 
different, ours being in general lower compared to the ones in |23j - this could be explained 
by the fact that these peaks appear at energy values that are well beyond the range of the 
validity of the approximation acj <^ 1 made in Ref. |2^; (iii) finally, differences appear 
also in the a* = case, with the emission curve dying out much faster in |23] than in our 
case. 



4.3 Angular momentum spectra 

We now turn to the rate of loss of the angular momentum of the higher-dimensional, 
rotating black hole through the emission of scalar fields on the brane. In Figs. ITUT a-b). 
this rate is shown, again, for the indicative cases of = 1 and n = 4, respectively, and 
for various values of the angular momentum parameter a^,. As is clear from both figures, 
the emission of angular momentum is clearly enhanced at all frequencies, as a* increases, 
and for all values of n. For low values of n, oscillations are again present but, as in the 
case of flux and power emission, there is an 'envelope' of Gaussian shape whose height 
increases and which peaks at higher values of uvh, as a=K increases. For higher values of 
n, the oscillations are smoothed out, and the enhancement is much more significant as 
a^, increases, due to the substantial increase both in the height and width of the emission 
curve. 




Figure 10: Angular momentum spectra for scalar particles from a rotating black hole, for 
(a) n = 1, and (b) n = 4, and various values of a*. 




Figure 11: Angular momentum spectra for scalar emission from a rotating black hole, for 
a* = 1 and various values of n. 



The dependence of the angular momentum loss rate on the number of additional 
spacelike dimensions n is more clearly depicted in Fig. ^2 for ^ fixed value of the angular 
momentum parameter, i.e a* = 1. As in the case of fiux and power spectra, a strong 
enhancement can be observed in the rate of loss of angular momentum of the black hole, 
as the dimensionality of spacetime increases. 

4.4 Total emissivities 

We finally turn to the computation of the total emissivities of particles, energy and an- 
gular momentum per unit time emitted by the black hole on the brane. These follow by 
integrating the quantities given in Eqs. (P3j) -(|7f j) with respect to ujTh, up to the value 
ujTh = 3. In Figs. IT^ a) and (b), we present two histograms depicting the dependence 
of the total fiuxes on the angular momentum parameter a^, (for fixed n = 4) and on the 
dimensionality of spacetime n (for fixed = 1), respectively. 

As a* increases, the histogram in Fig. IT^ a) shows that all fiuxes - particle, energy 
and angular momentum - are enhanced. More specifically, for n = 4, as a* goes from 
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Figure 12: Total emissivities for scalar emission on the brane from a rotating black hole 
function of (a) a*, for n = 4, and (b) n for a=K = 1. 



zero to 1.25, the number of particles emitted on the brane per unit time is enhanced 
by a factor of 2, the amount of energy emitted per unit time by a factor of 3, while, as 
expected, the angular momentum flux (which is zero when a=K = 0) increases significantly. 
The same enhancement appears in all three fluxes for all values of n, nevertheless the 
corresponding numbers are strongly n-dependent: for the case of n = 1, for instance, as 
a* goes again from zero to 1.25, the number of particles emitted per unit time increases 
by approximately 30%, the energy flux by a factor of 4, while the angular momentum 
flux, although significantly enhanced, is an order of magnitude smaller than the one for 
the n = 4 case. 

The dependence of the total emissivities on the dimensionality of spacetime is shown 
in the second histogram shown in Fig. ITW b). that clearly illustrates the enhancement of 
all fluxes as n increases. For a* fixed (a* = 1) and n varying from 1 to 7, the number of 
particles emitted by the black hole on the brane per unit time is enhanced by a factor of 
30, the energy flux by a factor of 100, and the angular momentum flux by a factor of 40. 

A note should be made at this point: the results presented above are only an ap- 
proximation to the exact, total emissivities that would follow by integrating the various 
fluxes over the complete frequency regime. In this work, numerical results have been 
produced only up to the frequency value of ur^ = 3 as the integration for higher values 
of uoTh takes an unrealistically long time. For small values of n and a^=, the flux dies away 
sufficiently quickly, therefore the contribution from uTh > 3 is negligible and the derived 
emissivities are an adequate approximation to the actual ones. However, as the value of 
either a^, or n increases, the peak of the Gaussian curve moves to higher energies, and the 
contribution of the part of the spectrum that we are missing to the exact emissivities is 
increasing. Nevertheless, the graphs presented in this subsection accurately describe the 
correct qualitative behaviour of the exact emissivities: the part of the spectrum that has 
been left out of our calculation would only increase further the enhancement of all three 
fluxes if it could be taken into account. 

We have also computed the integral over the frequency of the angular distribution 
of the particle flux (j28|) and power spectrum (j29|) to derive the emissivity solely as a 
function of the azimuthal angle 9. In Figs. IT^ a.b). we depict the power emissivity as a 
function of cos^ for fixed n and variable a*, and vice versa - the fiux emissivity exhibits 
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Figure 13: Power emissivity for scalar emission on the brane from a rotating black hole 
as a function of cos 6', (a) for n = 4 and variable a*, and (b) for a* = 1 and variable n. 



a similar behaviour and therefore is not shown here. The behaviour found agrees with 
the one deduced from the 3D graphs of subsections 4.1 and 4.2. For fixed n, as the 
angular momentum parameter a^, increases, the energy (and particle) flux is altogether 
enhanced and becomes more concentrated around the equatorial region; for fixed a*, the 
fluxes strongly increase as n increases, and so does again the proportion of the emission 
concentrated in the equatorial region. 

5 Discussion and Conclusions 

In this work we have performed a comprehensive analysis of the emission of Hawking 
radiation in the form of scalar fields from a (4 + n)-dimensional, rotating black hole on 
the brane. This analysis has provided complete results for the flux, energy and angular 
momentum spectra for the spin-down phase of the life of a higher-dimensional black hole, 
thus, filling a long-standing gap in the literature. By means of numerical analysis, we have 
been able to solve both the radial and angular equations for the scalar field modes, and 
our method is applicable for arbitrary energy of the emitted particles, angular momentum 
of the black hole, and number of extra, spacelike dimensions. 

We have first studied the fluxes of scalar particles, energy and angular momentum 
- integrated over the azimuthal angle 6 - for a range of values of the angular momentum 
of the black hole and number of extra dimensions n, in each function of the 

mode frequency uj. Due to their relation, a number of common features arise between the 
flux and power spectra: for low values of n, both emission rates are suppressed at the 
low energy regime but strongly enhanced in the intermediate- and high-energy one, as the 
angular momentum parameter increases; for high values of n, the enhancement is present 
at all frequency regimes. On the other hand, the rate of loss of the angular momentum 
of the black hole increases uniformly (i.e. at all frequency regimes) as a^, increases, for 
all values of n. When the angular momentum parameter is kept fixed and n increases, a 
clear, strong enhancement can be seen in all three spectra and in all frequency regimes. 
This enhancement can be easily explained by the fact that, for fixed a^,, as n increases, the 
temperature of the black hole, given by Eq. (fTTj) . increases too, thus leading to greater 
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emission rates. However, for fixed n, increasing a* lowers the temperature, and, in this 
case, the increased emission rates is due to enhanced super-radiance, which speeds up the 
spin- down process. 

We have also considered the angular distribution of the particle and energy fluxes. 
Again, a number of common features arise in the two spectra. An important feature is 
the dominance of the spherically-symmetric emission at the low-energy regime that leads 
to an almost ^-independent distribution even for a non-vanishing angular momentum 
parameter; as the energy of the emitted particles increases, however, a clear angular 
dependence - that is, the concentration of the emitted modes near the equatorial plane 
- appears in both spectra. We expect the angular variation in the emission rates to be a 
unique feature of the spectra during the spin-down phase. 

All spectra derived above, either integrated or not over the ^-coordinate, depend 
strongly on the energy regime that we look at, the value of the angular momentum 
parameter and the dimensionality of spacetime. In many cases, the various emission rates 
change considerably as the above parameters take different, low or high, values. In this 
work, care has been taken so that all relevant quantities, i.e. the angular eigenvalue 
as well as the radial and angular part of the scalar field, are computed via numerical 
analysis that is free from any assumptions and approximations that limit the values of 
the aforementioned three, fundamental parameters of the problem. 

Any assumptions made in this analysis involved only the mass of the produced 
black hole, being much larger than the fundamental Planck mass, and the size of the 
horizon of the black hole, being much smaller than the size of the extra dimensions. 
However, the values of and L themselves were never specified; they are free input 
parameters of the problem, and thus our analysis and corresponding results are valid 
for any gravitational theory with arbitrary number of extra dimensions and fundamental 
scale. Finally, although our analysis assumed that the extra spacetime is empty and thus 
fiat, our results are also valid in the case of a Randall- Sundrum ^U] type of black hole in 
the limit of a horizon radius being much smaller than the AdS radius. 

A number of interesting questions still remain open. Firstly, the calculation of the 
rate of spin-down of the black hole , and thus of the duration of this phase compared to 
the Schwarzschild one, requires, in addition to the brane emission, the scalar field radiation 
into the bulk for all values of the number of extra dimensions. The same calculation for 
the flux and energy spectra will reveal the amount of energy lost in the bulk and thus the 
remaining amount available for emission on the brane. Also, the remaining brane emission 
'channels', i.e. the ones for fermions and gauge bosons, need to be investigated too. It 
is thus necessary, and of great interest indeed, to extend our analysis for the emission of 
Hawking radiation from a rotating, higher-dimensional black hole for higher-spin particles. 
We plan to return to these questions in the near future. 
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